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Currently, molecular tunnel junctions are recognized as important active elements of various 
nanodevices. This gives a strong motivation to study physical mechanisms controlling electron 
transport through molecules. Electron motion through a molecular bridge is always somewhat 
affected by the environment, and the interactions with the invironment could change the energy 
of the traveling electron. Under certain conditions these inelastic effects may significantly modify 
electron transport characteristics. In the present work we describe inelastic and dissipative effects 
in the electron transport occurring due to the molecular bridge vibrations and stochastic thermally 
activated ion motions. We intentionally use simple models and computational techniques to keep 
a reader focused on the physics of inelastic electron transport in molecular tunnel junctions. We 
consider electron-vibron interactions and their manifestations in the inelastic tunneling spectra, 
polaronic effects and dissipative electron transport. Also, we briefly discuss long-range electron 
transfer reactions in macromolecules and their relation to the electron transport through molecular 
junctions. 
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I. Introduction 

Molecular electronics is known to be one of the most 
promising developments in nanoclcctronics, and the past 
decade has seen an extraordinary progress in this field 
(Aviram 2002, Cuniberti 2005, Nitzan 2001). Present 



activities on molecular electronics reflect the convergence 
of two trends in the fabrication of nanodevices, namely, 
the top-down device miniaturization through the litho- 
graphic methods and bottom-up device manufacturing 
through atom-engineering and self-assembly approach. 
The key element and basic building block of molecular 
electronics is a junction including two electrodes (leads) 
linked by a molecule, as schematically shown in the Fig. 
1. Usually, the electrodes are microscopic large but 
macroscopic small contacts which could be connected to 
a battery to provide the bias voltage across the junction. 
Such a junction may be treated as a quantum dot coupled 
to the charge reservoirs. The discrete character of energy 
levels on the dot (molecule) is combined with nearly con- 
tinuous energy spectra on the reservoirs (leads) occurring 
due to their comparatively large size. 

When the voltage is applied, an electric current flows 
through the junction. Successful transport experiments 
with molecular junctions (Ho 2002, Lortscher 2007, Park 
2000, Poot 2006, Reicher 2002, Smit 2002, Yu 2004) 
confirm their significance as active elements in nanode- 
vices. These include applications as rectifiers (molec- 
ular diodes), field effect transistors (molecular triods), 
switches, memory elements and sensors. Also, these ex- 
periments emphasize the importance of thorough analy- 
sis of the physics underlying electron transport through 
molecular junctions. Detailed understanding of the elec- 
tron transport at the molecular scale is a key step to 
future device operations. Theory of electron transport 
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FIG. 1: (Color online) Top panel: Schematic drawing of a 
junction including two electrodes and a molecule in between. 
Bottom panel: When the voltage is applied across the junc- 
tion electrochemical potentials and fiR differ, and the 
conduction window opens up. 



in molecular junctions is being developed in the last two 
decades, and main transport mechanisms are currently 
elucidated in general terms (Datta 2005, Imry 1999). 
However, progress of experimental capabilities in the field 
of molecular electronics brings new theoretical challenges 
causing further development of the theory. 

Speaking of transport mechanisms, it is useful to make 
a distinction between the elastic electron transport when 
the electron energy remains the same as it travels through 
the junction, and inelastic transport processes when the 
electron undergoes energy changes due to its interactions 
with the environment. There are several kinds of pro- 
cesses bringing inelasticity in the electron transport in 
mesoscopic systems including molecular junctions. Chief 
among these are electron-electron and electron-phonon 
scattering processes. These processes may bring signif- 
icant inelastic effects modifying transport properties of 
molecular devices and charging, desorption and chemical 
reactions as well. To keep this Chapter at a reasonable 
length we concentrate on the inelastic effects originating 
from electron-phonon interactions. 

In practical molecular junctions the electron transport 
is always accompanied by nuclear motions in the environ- 
ment. Therefore the conduction process is influenced by 
the coupling between electronic and vibrational degrees 
of freedom. Nuclear motions underlie the interplay be- 
tween the coherent electron tunneling through the junc- 
tion and inelastic thermally assisted hopping transport 
(Nitzan 2001). Also, electron-phonon interactions may 



result in polaronic conduction (Galperin 2005, Gutierrez 
2005, Kubatkin 2003, Ryndyk 2008), and they are di- 
rectly related to the junction heating (Segal 2003) and 
to some specific effects such as alterations in both shape 
of the molecule and its position with respect to the leads 
(Komeda 2002, Mitra 2004, Stipe 1999). The effects of 
electron-phonon interactions may be manifested in the 
inelastic tunneling spectrum which presents the second 
derivative of the current in the junction cPl/dV ver- 
sus the applied voltage V. The inelastic electron tun- 
neling microscopy has proven to be a valuable method 
for identification of molecular species within the con- 
duction region, especially when employed in combination 
with scanning tunneling microscopy and/or spectroscopy 
(Galperin 2004). 

Inelasticity in the electron transport through molec- 
ular junctions is closely related to the dephasing ef- 
fects. One may say that incoherent electron trans- 
port always includes an inelastic contribution with 
the possible exception of the low temperature range. 
The general approach to theoretically analyze electron 
transport through molecular junctions in the presence 
of dissipative/phase-breaking processes in both elec- 
tronic and nuclear degrees of freedom is based on 
the advanced formalisms ((Segal 2000, Skourtis 1995, 
Wingreen 1989,1993). These microscopic computational 
approaches have the advantages of being capable of pro- 
viding the detailed dynamics information. However, 
while considering stationary electron transport through 
molecular junctions, one may turn to less time consum- 
ing approach based on the scattering matrix formalism 
(Buttikcr 1986, Li 2001), as discussed below. 

II. Coherent transport 

To better show the effects of dissipation/dephasing 
on the electron transport through molecular junctions 
it seems reasonable to start from the case where these 
effects do not occur. So, we consider a molecule (pre- 
sented as a set of energy levels) placed in between two 
electrodes with nearly continuous energy spectra. While 
there is no bias voltage applied across the junction, the 
latter remains in equilibrium characterized by the equi- 
librium Fermi energy Ep, and there is no current flowing 
through it. When the bias voltage is applied, it keeps the 
left and right electrodes at different electrochemical po- 
tentials /il and /LtR. Then the electric current appears 
in the junction, and molecular energy levels located in 
between the electrochemical potentials fj,L and /ir take 
the major part in maintaining this current. Electrons 
from occupied molecular states tunnel to the electrodes 
in accordance with the voltage polarity, and the electrons 
from one electrode travel to another one using unoccu- 
pied molecular levels as intermediate states for tunneling. 
Usually, the electron transport in molecular junctions oc- 
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curres via highest occupied (HOMO) and lowest unoccu- 
pied (LUMO) molecular orbitals which work as channels 
for electron transmission. Obviously, the current through 
the junction depends on quality of contacts between the 
leads and the molecule ends. However, there also ex- 
ist the limit for the conductance in the channels. As 
was theoretically shown (Landauer 1970), the maximum 
conductance of a channel with a single spin degenerate 
energy level equals: 



G = — = (25. 8/cft)- 1 



(1) 



where e is the electron charge and h is the Planck's con- 
stant. This is a truly remarkable result for it proves that 
the minimum resistance i?o = Gq 1 of a molecular junc- 
tion cannot become zero. In another words one never can 
short-circuit a device operating with quantum channels. 
Also, the expression (JIJ shows that the conductance is a 
quantized quantity. 

Conductance g in practical quantum channels associ- 
ated with molecular orbitals can take on values signif- 
icantly smaller that Go, depending on the derealiza- 
tion in the molecular orbitals participating in the elec- 
tron transport. In molecular junctions it also strongly 
depends on the molecule coupling to the leads (quality 
of contacts), as was remarked before. The total resis- 
tance r — g^ 1 includes contributions from contact and 
molecular resistances, and could be written as (Wingreen 
1993): 



Gn 



1 



1 -T 
T 



(2) 



Here, T is the electron transmission coefficient which 
generally takes on values less than unity. 

The general expression for the electric current flowing 
through the molecular junction could be obtained if one 
calculates the total probability for an electron to travel 
between two electrodes at a certain tunnel energy E and 
then integrates the latter over the whole energy range 
(Datta 1995). This results in the well known Landauer 
expression: 



I = 



Tit) 



T(E)[f L (E)-f R (E)]dE. 



(3) 



Here, fL,R.(E) are Fermi distribution functions for the 
electrodes with chemical potentials hl,r, respectively. 
The values of /j.l,r differ from the equilibrium energy 
Ep , and they are determined with the voltage distri- 
bution inside the system. Assuming that the coherent 
tunneling predominates in the electron transport, the 
electron transmission function is given by (Datta 1997, 
Samanta 1996): 



T{E) = 2Tr{A L GA R G + } 



(4) 



where the matrices A^.^ represent imaginary parts of 
self-energy terms ^l,r describing the coupling of the 



molecule to the electrodes, and G is the Green's function 
matrix for the molecule whose matrix elements between 
the molecular states (z| and \j) have the form: 



Gi 



(i\E-H\j). 



(5) 



Here, H is the molecular Hamiltonian including self- 
energy parts S L ,fl. 

When a molecule contacts the surface of electrodes, 
this results in charge transfer between the molecule and 
the electrodes, and in modification of the molecule en- 
ergy states due to redistribution of the electrostatic po- 
tential within the molecule. Besides, the external voltage 
applied across the junction brings additional changes to 
the electrostatic potential further modifying molecular 
orbitals. The coupling of the molecule to the leads may 
also depend on the voltage distribution. So, generally, 
the electron transmission T inserted in the Eq. ([3]) and 
electrochemical potentials (J>l,r depend on the electro- 
static potential distribution in the system. To find the 
correct distribution of the electric field inside the junction 
one must simultaneously solve the Schrodinger equation 
for the molecule and the Poisson equation for the charge 
density, following a self-consistent converging procedure. 
This is a nontrivial and complicated task, and signifi- 
cant effort was applied to study the effect of electrostatic 
potential distribution on the electron transport through 
molecules (Damle 2001, Di Ventra 2000, Galperin 2006, 
Lang 2000, Mujica 2000, Xue 2001,2003,2004). Here we 
put these detailed considerations aside, and we use the 
simplified expression for the electrochemical potentials: 



H L = E F + r?|e|V; 



E f 



(l-v)\e\V (6) 



where the parameter rj indicates how the bias voltage 
is distributed between the electrodes. Also, we assume 
that inside the molecule the external electrostatic field is 
screened due to the charge redistribution, and the elec- 
tron transmission is not sensitive to the changes in the 
voltage V. Although very simple, this model allows to 
analyze the main characteristics of the electron transport 
through molecular junctions. Within this model one may 
write down the following expression for the self-energy 
parts ((D'Amato 1990): 



0)ij 



E 



E 



- £k,0 + 1 



(7) 



Here, /3 G L, R, Tik,/3 is the coupling strength between 
"i" -th molecular state and "fc" -th state on the left / right 
lead, efe :i s are the energy levels on the electrodes, and s 
is a positive infinitesimal parameter. Assuming that the 
molecule is reduced to a single orbital with the energy 
Eo (a single-site bridge), the Green's function accepts 
the form: 



G(E) = 



1 



E - E — E L — Sj 



(8) 
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FIG. 2: (Color online) Coherent electron transmission (left 
panel) and current (right panel) versus bias voltage applied 
across a molecular junction where the molecule is simulated 
by a single electronic state. The curves are plotted assuming 
Al = Ajj = O.leV, E = -0.5eV, T = 30K. 



FIG. 3: (Color online) Coherent electron transmission (left 
panel) and current (right panel) through a junction with a 
five electronic states bridge. The curves are plotted for Al = 
Ar = O.leV, b = O.SeV, E = -0.5eV, T = 30K. 



Accordingly, in this case one may simplify the expression 
(UJ for the electron transmission: 



T(E) 



(E - E y + (A L + A R y 



(9) 



To elucidate the main features of electron transport 
through molecular junctions we consider a few exam- 
ples. In the first example we mimic a molecule as a one- 
dimensional chain consisting of N identical hydrogen- 
like atoms with nearest neighbors interaction. We as- 
sume that there is one state per isolated atom with the 
energy E , and that the coupling between the neighbor- 
ing sites in the chain is characterized by the parameter 
b. Such a model was theoretically analyzed by D'Amato 
and Pastawski and in some other works (see e.g. Mujica 
1994). Basing on the Eqs. Q, (P, it could be shown that 
for a single-site chain (N = 1) the electron transmission 
reveals a well distinguished peak at E = Eq, , shown in 
the Fig. 2 The height of this peak is determined with the 
coupling of the bridge site to the electrodes. The peak 
in the electron transmission arises because the molecular 
orbital E = E works as the channel/bridge for electron 
transport between the leads. Similar peaks appear in the 
conductance g = dl/dV. Assuming the symmetric volt- 
age distribution (rj = 1/2) the peak in the conductance 
is located at V = ±2E$. As for the current voltage char- 
acteristics, they display step-like shapes with the steps at 
V = ±2Eq. When the chain includes several sites we ob- 
tain a set of states (orbitals) for our bridge instead of the 
single state E = Eq, and their number equals the num- 
ber of sites in the chain. All these states are the channels 
for the electron transport. Correspondingly, the trans- 
mission reveals a set of peaks as presented in the Fig. 3. 
The peaks are located within the energy range with the 
width 46 around E = Eq. The coupling of the chain 
ends to the electrodes affects the transmission, especially 



near E = Eq. As the coupling strengthens, the transmis- 
sion minimum values increase. Now, the current voltage 
curves exhibit a sequence of steps. The longer is the chain 
the more energy levels it possesses, and the more steps 
appear in the I — V curves. 

The second example concerns the electron transport 
through a carbon chain placed between copper elec- 
trodes. In this well as for practical molecules a 
preliminary step in transport calculations is to compute 
the relevant molecular energy levels and wave functions. 
Usually, these computations are carried out employing 
quantum chemistry software packages (e.g. GAUSSIAN) 
or density-functional based software. Also, a proper 
treatment of the molecular coupling to the electrodes is 
necessary for it brings changes into the molecular energy 
states. For this purpose one may use the concept of an 
"extended molecule" proposed by Xue, Datta and Rat- 
ner (Xue 2001). The point of this concept is that only 
a few atoms on the surface of the metallic electrode are 
significantly disturbed when the molecule is attached to 
the latter. These atoms are located in the immediate 
vicinity of the molecule end. Therefore one may form 
a system consisting of the molecule itself and the atoms 
from the electrode surfaces perturbed by the molecular 
presence. This system is called the extended molecule 
and treated as such while computing the molecular or- 
bitals. In the considered example the extended molecule 
included four copper atoms at the each side of the car- 
bon chain. The results for electron transport are shown 
in the Fig. 4. Again, we observe a comb-like shape of the 
electron transmission corresponding to the set of trans- 
port channels provided by the molecular orbitals, and the 
stepwise I — V curve originating from the latter. Trans- 
port calculations similar to those described above were 
repeatedly carried out in the last two decades for vari- 
ous practical molecules (see e.g. Xue 2001,2003, Galperin 
2006, Zimbovskaya 2002,2003,2008). 
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FIG. 4: (Color online) Coherent electron transmission (left 
panel) and current (right panel) through a carbon chain cou- 
pled to the copper leads at T = 30K. 



FIG. 5: (Color online) Schematic drawing illustrating inelas- 
tic electron transport through a molecular junction within the 
Buttiker model. 



III. Buttiker model for inelastic transport 

An important advantage of the phenomenological 
model for the incoherent/inelastic quantum transport 
proposed by Buttiker (Buttiker 1986) is that this model 
could be easily adapted to analyze various inelastic effects 
in the electron transport through molecules (and some 
other mesoscopic systems) avoiding complicated and time 
consuming advanced methods such as based on the non- 
equilibrium Green's functions formalism (NEGF). 

Here, we present the Buttiker model for a simple 
junction including two electrodes linked by a single-site 
molecular bridge. The bridge is attached to a phase- 
randomizing electron reservoir, as shown in the Fig. 5. 
Electrons tunnel from the electrodes to the bridge and 
vice versa via the channels 1 and 2 . While on the 
bridge, an electron could be scattered into the channels 3 
and 4 with a certain probability e. Such electron arrives 
at the reservoir where it undergoes inelastic scattering 
accompanied by phase-breaking and then the reservoir 
reemits it back to the channels 3 and 4 with the same 
probability. So, within the Buttiker model the electron 
transport through the junction is treated as the combi- 
nation of tunnelings through the barriers separating the 
molecule from the electrodes and interaction with the 
phase-breaking electron reservoir coupled to the bridge 
site. The key parameter of the model is the probability 
e which is closely related to the coupling strength be- 
tween the bridge site and the reservoir. When e = the 
reservoir is detached from the bridge, and the electron 
transport is completely coherent and elastic. Within the 
opposite limit (e = 1) electrons are certainly scattered 
into the reservoir which results in the overall phase ran- 
domization and inelastic transport. 

Within the Buttiker model the particle fluxes outgoing 
from the junctions J- could be presented as linear com- 
binations of the incoming fluxes Jk where the indexes 



i, k label the channels for the transport: 1 < i, k < 4. 



J'i — ^ TikJk- 



(10) 



The coefficients in these linear combinations are ma- 
trix elements of the transmission matrix which are re- 
lated to the elements of the scattering matrix S, namely: 
Tik = \Sik\ 2 - The matrix S expresses outgoing wave 
amplitudes b[ ,b' 2l a' 3 , a' A in terms of the incident ones 
^i7&2, 13,04- To provide the charge conservation in the 
system, the net current in the channels 3 and 4 linking 
the system with the dephasing reservoir must be zero, so 
we may write: 



J3 + J4- 4 -J' 4 = 0. 



(11) 



The transmission for the quantum transport could be 
defined as the ratio of the particle flux outgoing from the 
system and that one incoming to the latter. Solving Eqs. 
(jTO]) . (fTTj) we obtain: 



T{E) = 3 1 



T 



21 



K X -K 2 

2R 



(12) 



where: 



Ki = T31 + Tn] K 2 = T 23 + T 24 ; 
R = I33 + + T43 + T34 . 



(13) 



For the junction including the single-site bridge the scat- 
tering matrix S has the form (Buttiker 1986): 



In 



a 2 r 2 



0*1*2 
\ afitir-i 



(3t x aptir 2 \ 
afirit 2 j3t 2 
ft 2 ri arir 2 — a 
OLT\r 2 — a P 2 r 2 J 
(14) 

Here, Z = 1 — a 2 rir 2 , a = VI — e , P = \[t and ri 2 
and ii 2 are the amplitude reflection and transmission 
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FIG. 6: (Color online) Electron transmission (left panel) and 
current (right panel) computed within the Buttiker model at 
various values of the dephasing parameter e, namely: e = 
(dotted lines), e — 0.5 (dashed lines), and e = 1 (solid 
lines). The curves are plotted assuming that the molecule is 
simulated by a single orbital with Eo = — 0.5eV at Al ~ 
A R = Q.leV, T = 30K. 



coefficients for the two barriers. Later the expression for 
this matrix suitable for the case of multi-site bridges in- 
cluding several inelastic scatterers was derived (Li 2001). 

Assuming for certainty the charge flow from the left 
to the right we may write down the following expression 
(Zimbovskaya 2005): 

T(E) = g(E)(l + a 2 )[ 9 (E)(l + a 2 ) + l~a 2 ] 
[g(E)(l - a 2 ) + 1 + a 2 ] 2 

where 

M- ^B-w^ Tzii?- (16) 

Now, the electron transmission strongly depends on the 
dephasing strength e. As shown in the Fig. 6 coherent 
transmission (e = 0) exhibits a sharp peak at E = Eq 
which gives a step-like shape to the volt-ampere curve, 
as was discussed in the previous section. In the presence 
of dephasing the peak gets eroded. When the e value 
approaches 1 the I — V curve becomes linear, corrobo- 
rating the ohmic law for the inelastic transport. 

Within the Buttikcr's model e is introduced as a phe- 
nomenological parameter whose relation to the micro- 
scopic characteristics of the dissipative processes affect- 
ing electron transport through molecular junctions re- 
mains uncertain. To further advance this model one 
should find out how to express e in terms of the relevant 
microscopic characteristics for various transport mecha- 
nisms. This should open the way to make a link between 
the phenomenological Buttiker model and NEGF. Such 
attempt was carried out in recent works (Zimbovskaya 
2005,2008) where the effect of stochastic nuclear motion 
on the electron transport through molecules was ana- 
lyzed. 



IV. Vibration-induced inelastic effects 

Interaction of electrons with molecular vibrations is 
known to be an important source of inelastic contribution 
to the electron transport through molecules. Theoret- 
ical studies of vibrationally inelastic electron transport 
through molecules and other similar nanosystems (e.g. 
carbon nanotubes) were carried out over the past few 
years by a large number of authors (Cornaglia 2004,2005, 
Donarini 2006, Egger 2008, Galperm 2007, Gutierrez 

2006, Kushmerick 2004, Mii 2003, Rynduk 2007, Siddiqui 

2007, Tikhodeev 2004, Troisi 2006, Zazunov 2006, 2007, 
Zimmerman 2008). Also, manifestations of the electron- 
vibron interactions were experimentally observed (Agrait 

2003, Djukic 2005, Lorente 2001, Qiu 2004, Repp 2005, 
Segal 2001, Smit 2004, Tsutsui 2006, Wang 2004, Wu 

2004, Zhitenev 2002). To analyze vibration induced ef- 
fects on the electron transport through molecular bridges, 
one must assume that molecular orbitals are coupled to 
the phonons describing vibrations. While on the bridge, 
electrons may participate in the events generated by their 
interactions with vibrational phonons. These events in- 
volve virtual phonon emission and absorption. For rather 
strong electron-phonon interaction this leads to the ap- 
pearance of metastable electron levels which could par- 
ticipate in the electron transport through the junctions, 
bringing an inelastic component to the current. As a re- 
sult, vibration induced features occur in the differential 
molecular conductance dl/dV and in the inelastic tun- 
neling spectrum d 2 I /dV 2 . This was observed in the ex- 
periments (see e.g. Qiu 2004, Zhitenev 2002). Sometimes 
even current voltage curves themselves exhibit extra step 
originating from the electron- vibron interactions (Djukic 
2005). 

Particular manifestations of electron-vibrionic effects 
in the transport characteristics are determined by the 
relation of three relevant energies. These are the cou- 
pling strengths of the molecule to the electrodes &-l,r, 
the electron-phonon coupling strength A and the ther- 
mal energy kT (k is the Boltzman constant). When the 
molecule is weakly coupled to the electrodes (A l ,r <c 
A, Ml) and the temperature is low (kT <C A L ^), the 
electron transfer through the junction may give rise to a 
strong vibrational excitation, and one may expect pro- 
nounced vibrational resonance structure in the electron 
transmission to appear. Correspondingly, extra steps 
should occur in the I— V curves. Proper theoretical con- 
sideration of the electron transport in this regime could 
be carried out within the approach proposed by Wingreen 
et al (Wingreen 1989). Here, we employ a very simple 
semi-quantitative approximation which, nevertheless, al- 
lows to qualitatively describe this structure including the 
effect of higher phonon harmonics on the transport char- 
acteristics. 

We consider a junction including a single-site bridge 
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FIG. 7: (Color online) The electron transmission versus en- 
ergy (left panel) and the conductance versus voltage (right 
panel) for a junction with the molecular bridge simulated 
by a single electronic state weakly coupled to the leads: 
A L = A R = O.OleV; E = 0.5eV. Solid lines are plotted 
assuming that the bridge is coupled to a single phonon mode 
(Ml — 0.22eV, A = 0.3eV). Dashed lines correspond to the 
coherent electron transport. 
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FIG. 8: (Color online) The inelastic electron tunneling spec- 
trum plotted against the bias voltage at HQ, = 0.22eV, A = 
0.3eV. Left panel: Al = Ar = O.OleV, E = 0.5eV. 
Right panel: A L = A R = 0.5eV, E = 0.5eV (solid line), 
E Q = ~0.5eV (dashed line). 



which is coupled to a single vibrational mode with the 
frequency VL. An electron on the bridge may virtually 
absorb several phonons which results in the creation of 
a set of metastable states with the energies E n = Eq + 
nMl (n = 0,1,2, ...). Here, the energy Eq is shifted with 
respect to Eq due to the electron-phonon interaction. 
The difference in these energies E p is called a polaronic 
shift and could be estimated as E p — A 2 /Ml (Gutierrez 
2006, Ryndyk 2007, Wingreen 1989). At weak coupling 
of the bridge state to the electrodes the lifetime of these 
metastable states is long enough for them to serve as 
channels for the electron transmission. Therefore, one 
may roughly approximate the transmission as a sum of 
contributions from all these channels. The terms in the 
sum have the form similar to the well-known expression 
for the coherent transmission (see Eq. (j9])). However, 
every term includes the factor P(n) which corresponds 
to the probability of the metastable state to appear. So, 
we obtain: 



T{E) = 2A L A R J2 



P(n) 



Here (Cizek 2004): 



[E — Eq — nMiy + (A L + A R ) 2 

(17) 



1 ( X% 



(18) 



The phonon induced peaks in the transmission are dis- 
played in the Fig. 7 along with the transmission peak for 
the coherent transport through a single-site bridge. As 
expected, the coupling of electronic degrees of freedom to 
the vibrational motion splits the single peak in the coher- 
ent transmission into the set of smaller peaks associated 
with vibrational levels. The peaks could be resolved when 



A_l,j?. <C Ml. This agrees with the results of the earlier 
theoretical works (Wingreen 1989) as well as with the ex- 
periments (Qiu 2004, Zhitenev 2002). Phonon-induced 
peaks in the transmission give rise to the steps in the 
I —V curves and rather sharp features (peaks and dips) 
in the inelastic tunneling spectrum (IETS). The latter 
are shown in the Fig. 8 (left panel) and they resemble 
those obtained using proper NEGF based calculations 
(Galperin 2004). 

When the molecule is strongly coupled to the elec- 
trodes {Al.r ^S> A) and the temperature is still low 
(kT <C A, Ml), electron-vibrionic interaction effects 
are less pronounced. Both current and conductance 
are weakly affected by the electron-phonon coupling 
(Galperin 2007, Tikhodeev 2004). However, the IETS 
features remain distinguishable. These features appear 
at the threshold V — Ml/\e\ which corresponds to the 
opening of a channel for inelastic transport. To ana- 
lyze IETS in a simplest way one may use the result for 
electron transmission derived within the Buttiker model 
where the dephasing parameter e is expressed in terms 
of the relevant energies, namely: 



pi i 



2{A L + A R ) + T ph 



(19) 



where T p h = 2Im(E p /j), end is the self-energy term 
originating from the electron-phonon interaction. Basing 
on the nonequilibrium Green's function formalism, the 
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expression for T p h was derived in the form (Mii 2003): 

T ph (E) = 2nX 2 J du)p{u){N{u)\p el {E - hu>) + Pel (E + Hu)} 

+ (1 + N(w))([l - n(E - Huj)]p e i(E - hu)) 

+ n(E + hw)p e i(E + hui)) 

+ N(u)([l - n{E + Hu>)]p el {E + tuS) 

+ n(E -Hu)p e i(E -hu>))}. (20) 

Here, Pei(E) and p P h{u) are the phonon and elec- 
tron densities of states, respectively; N{ui) is the Bose- 
Einstein distribution function at the temperature T, 
and: 



i(E) = ±[f L (E) + f R (E)]. 



(21) 



While considering a junction with a single-site bridge 
state coupled to the single vibrational mode, we can write 
the following expressions for p p h and p e i : 



Pph{u) 



1 
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Pel(E) = 



7T {E - E ) 2 + (A L + A R ) 



(22) 



(23) 



where the polaron shift is neglected, and the parameter 
7 characterizes the broadening of the maximum in p p h 
at oj = Q due to the interaction of the vibrionic mode 
with the environment. At low temperatures we may sig- 
nificantly simplify the expression for r„/j. Within the 
conduction window pr < E < p^ we get: 

T p h{E) xnrX 2 ^ J dwp p h(<jj)p e i(E + tvjj) 




(E-fi R )/h 



dujp ph (uj)p e i(E - Tiki) > . (24) 



Omitting from consideration the coupling of the phonon 
mode to the environment (7 — ¥ 0) we may easily carry 
out integration over uj in the Eq. (1241) , and we arrive at 
the result: 



T ph {E) ^n\ 2 {p el {E + m)0(ji L -m-E) 
+ Pel (E - Huj)6(E -PR- m)} 



(25) 



where 9{x) is the step function. 

Substituting the approximation for T p h given by Eq. 
(f2"4")l or Eq. (|2l)|) into the expression (fTO)) for the dephas- 
ing parameter e and employing the earlier result (|15[) 
for the electron transmission, we may calculate transport 
characteristics. The adopted simplified approach, as well 
as NEGF based calculations, shows that the IETS signal 
appears at the threshold determined by the frequency of 
the vibrational mode (see Fig. 8, right panel). At first 



glance one may expect the net current through the junc- 
tion to increase at the threshold. Indeed, the inelastic 
contribution to the current increases from zero to a cer- 
tain nonzero value at this threshold for the channel for 
inelastic transport opens up. However, more thorough 
studies show that both elastic and nonelastic contribu- 
tions to the net current undergo changes at the inelas- 
tic tunneling threshold, and the elastic current could de- 
crease there, as was first shown by Persson and Baratoff 
(Persson 1987). Moreover, this decrease in the elastic 
current may overweigh the contribution coming from the 
inelastic channel. Depending on the relative value of the 
elastic and inelastic contributions to the net current near 
threshold, the IETS reveals a peak or a dip at the corre- 
sponding voltage. Experiments corroborate the variety in 
the IETS taken for molecular junctions (Hahn 2000, Dju- 
kic 2005, Wang 2004, Zhitenev 2002). The shape of the 
signal is very sensitive to the characteristics of the junc- 
tion such as the position of the electronic state, electron- 
phonon and molecule-to leads coupling strengths, and 
the vibrionic frequency. For instance, at a very strong 
coupling of the molecular bridge to the leads, it could so 
happen that the backscattered by the negatively biased 
electrode electrons whose energies belong to the conduc- 
tion window between p^ and and pn, are locally de- 
pleted near the junction. In this situation the opening of 
the inelastic channel may cause the increased reflection 
(otherwise forbidden by the Fermi exclusion) leading to 
the decrease in the conduction. Consequently, the de- 
scribed scenario should result in the dip in IETS signal. 
Also, as concluded in the recent work of Ryndyk and 
Cunibcrti (Ryndyk 2007) the above discussed sideband 
phonon-induced features in the electron spectral density 
could give rise to the corresponding features in the differ- 
ential conductance dl/dV and IETS assuming that the 
molecule coupling to the leads is not too strong. Con- 
tributions from these sideband features may be respon- 
sible for the shape of the IETS signal at the threshold 
of the inelastic tunneling channel. These contributions 
could produce an extra inelastic signal, as well. The lat- 
ter appears as an additional peak or dip in the differential 
conductance. 

The question of current decrease/increase at the 
phonon excitation threshold which corresponds to the 
peak/dip in the IETS is not completely answered so 
far, and the appropriate theory is being developed (Bal- 
seiro 2006, Egger 2008). Nevertheless, it is presently un- 
derstood that three relevant energies, namely, molecule- 
electrode couplings Al^r, electron-phonon coupling 
strength A and the phonon energy Ml play very im- 
portant part in determining the shape of the IETS sig- 
nal. Varying these parameters one may convert the IETS 
signal from a peak to a dip and vice versa. 

At finite temperatures molecular vibrations always oc- 
cur in the presence of stochastic nuclear motions. These 
motions could be described as a phonon thermal bath. 
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Coupling of vibrational modes to this bath further af- 
fects the electron transport causing energy dissipation. 
The dissipative processes must be taken into account to 
properly analyze the effects of electron-phonon interac- 
tions in the electron transport. Also, the displacements 
of ions involved in the molecular vibrations are accompa- 
nied with the changes in the electrostatic field inside the 
molecule. This could give rise to polaronic effects in the 
electron transport. We discuss these issues in the next 
sections. 



V. Dissipative transport 

Electron transport through molecular junctions is al- 
ways accompanied by stochastic nuclear/ion motions in 
the close environment. Interactions of travelling electron 
with these environmental fluctuations causes energy dis- 
sipation. The importance of dissipative effects depends 
on several factors. Among these factors the temperature 
and the size and complexity of the molecular bridge are 
predominating. The temperature determines the inten- 
sity of the nuclear motions, and the size of the molecule 
determines so-called contact time, that is the time for 
electron to travel through the junction and, in conse- 
quence, to contact the environment. It was shown by 
Buttiker and Landauer (Buttiker 1985) that the contact 
time is proportional to the number of sites (subunits) in 
the molecule providing intermediate states for the elec- 
tron tunneling. For small molecules at low temperatures 
the contact time is shorter than characteristic times for 
fluctuations in the environment, so, the effect of the latter 
on the electron transport is not very significant. One may 
expect small broadening of the molecule energy states to 
occur which brings a moderate erosion of the peaks in 
the electron transmission and steps in the I — V char- 
acteristics. 

On the contrary, in large sized molecules such as pro- 
teins and DNA, electron transport is accompanied by 
strong energy dissipation. The significance of the system- 
environment interactions in macromolecules was recog- 
nized long ago in studies of long range electron transfer 
reactions. In these reactions electrons travel between dis- 
tant sites on the molecule called a donor and an acceptor. 
It was established that when an electron initially local- 
ized on the donor site moves to the acceptor site with 
a lower energy, the energy difference must be dissipated 
to the environment to provide the irreversibility of the 
transfer (Garg 1985). 

A usual way to theoretically analyze dissipative effects 
in the electron transport through molecules is to intro- 
duce a phonon bath representing the random motions 
in the environment. In general, there is no one/several 
dominating modes in the bath. Instead, the bath is char- 
acterized by a continuous spectral function of all relevant 
phonon modes p p h(uj)- The electrons are supposed to be 



coupled to the phonon bath, and this coupling is spec- 
ified by the spectral function. The particular form of 
Pp/i(w) may be found basing on the molecular dynamics 
simulations. However, to qualitatively study the effect of 
dissipation on the electron transport one can employ the 
expression (Mahan 2000): 



P P h( 



X — exp 

Ul c \ LJ. 



(26) 



where the parameter A characterizes electron-phonon 
coupling strength, and ui c is the cutoff frequency for the 
bath related to its thermal relaxation time r c = w~ . 

To illustrate possible effects of dissipation on the elec- 
tron transport we return back to our simple model 
where the molecular bridge is represented by a single 
state. Now, we assume that this state is coupled to 
the phonon bath. This model is hardly appropriate 
to properly analyze dissipative effects in the electron 
transport through practical molecules for the molecule 
length is very important for dissipative effects to be 
pronounced. Nevertheless, it still could serve to basi- 
cally outline main features of the dissipative electron 
transport through molecular junctions. Also, the pro- 
posed model could be useful to analyze electron transport 
in doped polyacetylene/polyaniline-polyethylene oxide 
nanofibers (Zimbovskaya 2008). These conducting poly- 
mers could be treated as some kind of granular metals 
where metallic-like regions (grains) are embedded in the 
poorly conducting medium of disorderly arranged poly- 
mer chains (MacDiarmid 2001). While in metallic state, 
the intergrain electron transport in these nanofibers is 
mostly provided by electron tunneling through interme- 
diate states on the polymer chains between the grains 
(Prigodin 2002). In this case the contact time could be 
long enough for the effects of dissipation to be well man- 
ifested which justifies the adoption of the above model. 

Again, one may carry out transport calculations using 
Eq. (|I5p for the electron transmission, and expressing 
e in terms of of the relevant energies. Substituting Eq. 
(|2"6l in the expression (j2T)]) we may calculate T p h(E). The 
energy dissipation effects are more distinctly prononced 
at moderately high temperatures, so we assume kT 
Huj c . Then the main contribution to the integral over oj 
in the Eq. (|20j) originates from the low frequency region 
lj <C w c ; and we obtain the following approximation: 



T ph (E) 



8kTXT 



(E - £ )2 + T2 



(27) 



where T = + Ar + ^Y p h(E). Solving the obtained 
equation for Y p h and using the Eq. (| 19|) we obtain (Zim- 
bovskaya 2008): 



p 2 (i + v /iT7) 



E-E 
A L + A, 



5(1 + ^ 



(28) 
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FIG. 9: (Color online) Schematic drawing of a molecular junc- 
tion where the molecular bridge is coupled to the phonon bath 
via side chains. 



where p 2 = 32kT\/(A L + A R ) 2 . 

Voltage dependencies of the conductance computed 
basing on Eqs. (O, (H"5j) , (|2"51) are presented in the Fig. 
10 (left panel). One may see that at low values of the bias 
voltage the electrons coupling to the phonon bath brings 
an enhancement in the conduction. The effect becomes 
reversed as the voltage grows above a certain value. This 
happens because the phonon-induced broadening of the 
molecular level (the bridge) assists the electron trans- 
port at small bias voltage. When the voltage rises this 
effect is surpassed by the scattering effect of phonons 
which resists the electron transport. The significance of 
the electron-phonon interactions is determined by the ra- 
tio of the coupling constant A and the self-energy terms 
describing the bridge coupling to the electrodes. The 
phonon bath makes a significant effect on the transport 
characteristics when A > Al,r- 

A dissipative electron transport through large DNA 
molecules was studied both theoretically and experimen- 
tally (Gutierrez 2005, Xu 2004). Theoretical studies were 
based on a model where the molecule was simulated by a 
tight-binding chain of sites linking the electrodes, and at- 
tached side chains. Electrons are allowed to travel along 
the bridge chain and to hop to the near by side chains. 
These chains are coupled to the phonon bath providing 
the energy dissipation (see Fig. 9). Although proposed 
for the specific kind of poly (dG) -poly (dC) molecules, 
this model seems to be quite generic and useful for a 
larger class of macromolecules. 

Several coupling regimes to the bath may be analyzed. 
The most preferred regime for dissipative effects to ap- 
pear is the strong-coupling limit defined by the condition 
\/lj c > 1. Within this regime the characteristic time for 
the electron bath interactions is much shorter than typ- 
ical electron time scales. Consequently, the bath makes 
a significant impact on the molecule electronic struc- 
ture. New bath- induced states appear in the molecular 
spectrum inside the HOMO-LUMO gap. However, these 
states are strongly damped due to the dissipative action 
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FIG. 10: (Color online) Left panel: The electron conduc- 
tance versus voltage for a junction with a single electronic 
state bridge directly coupled to the phonon thermal bath. 
The curves are plotted assuming Eo — OAeV, Al — Ar = 
O.leV, T = 30K, A = 0.3eV (solid (line), A = 0.05eV 
(dashed line). Right panel: Electron transmission through 
the junction in the case when the bridge state interacts with 
the phonon bath via the side chain coupled to the bridge state 
with the coupling parameter w. The curves are plotted as- 
suming Al = Ar = A, it;/ A = 20, E = 0, A = 0.3eV 
(solid line), A = 0.05eV (dashed line). 



of the bath (Gutierrez 2005). As a result, a small finite 
density of phonon-induced states appears inside the gap 
supporting electron transport at low bias voltage. So 
again, the environment induces incoherent phonon as- 
sisted transport through molecular bridges. For illus- 
tration we show here the results of calculations carried 
out for a toy model with a single-site bridge with a side 
chain attached to the latter. The side chain is supposed 
to be coupled to the phonon bath. The results for the 
electron transmission are displayed in the Fig. 10 (right 
panel). We see that the original bridge state at E = is 
completely damped but two new phonon-induced states 
emerge nearby which could support electron transport. 

An important characteristics of the dissipative electron 
transport through molecular junctions is the power loss 
in the junction, that is the energy flux from electronic 
into phononic system. Assuming the current flow from 
the left to the right, this quantity may be estimated as 
the sum of the energy fluxes Ql,r at the left and right 
terminals (leads): 



P = Ql + Qr- 



(29) 



One may express the energy fluxes in terms of renormal- 
ized currents at the electrodes Il,r{E) which are defined 
as follows (Datta 2005): 



Il.r = I L , R (E)dE. 

Then Ql,r may be presented in the form: 
1 



!L,R 



EI L , R {E)dE. 



(30) 



(31) 
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We remark that the current Il,r in the Eq. (f30l) arc 
related to the corresponding leads, and their signes are 
accordingly defined. An outgoing current is supposed 
to be positive whereas an incoming one is negative for 
each lead. To provide electric charge conservation in the 
junction one must require that I = II = —Ir for the 
chosen direction of the current flow. Therefore the en- 
ergy fluxes also have different signs. As for the Ql,r 
magnitudes, they may differ only if the renormalized cur- 
rents Il(E) and Ir(E) are distributed over energies in 
different ways. This cannot happen in the case of elas- 
tic transport for in this case Il(E) = —Ir(E). How- 
ever, if the transport process is accompanied with the 
energy dissipation, the energy distributions for II and 
Ir may differ. In this case electrons lose some energy 
while moving through the junction and this gives rise to 
the differencies in the renormalized currents energy dis- 
tributions. For instance, in the case when the electrodes 
are linked with a single-state bridge, the energy distribu- 
tion of Il(E) has a single maximum whose position is 
determined by the site energy Eq and the applied bias 
voltage V. Assuming that the average energy loss due 
to dissipation could be estimated as AE, the maximum 
in the Ir(E) distribution is shifted by this quantity, so 
the current Ir(E) flows at lower energies compared to 
Il{E) . This results in the power loss and Joul heating 
in the junction (Segal 2003). 

VI. Polaron effects: Hysteresis, switching and 
negative differential resistance 

While studying electron transport through molecular 
junctions hysteresis in the current- voltage characteristics 
was reported in some systems (Li 2003). Multistability 
and stochastic switching were reported in single-molecule 
junctions (Lortscher 2007) and in single metal atoms cou- 
pled to a metallic substrate through a thin ionic insulat- 
ing film (Olsson 2007, Repp 2004). 

Coupling of an electron belonging to a certain atomic 
energy level to the displacements of ions in the film brings 
a possibility of polaron formation in there. This leads to 
the polaron shift in the electron energy. It was noticed 
that multistability and hysteresis in molecular junctions 
mostly occurred when the molecular bridges included 
centers of long-living charged electronic states (redox cen- 
ters). On these grounds it was suggested that hystere- 
sis in the I — V curves observed in molecular junctions 
appear due to formation of polarons on the molecules 
(Galperin 2005). 

The presence of the polaron shift in the energy of a 
charged (occupied) electron state creates a difference be- 
tween the latter and the energy of the same state while 
it remains unoccupied. Assuming for simplicity a single- 
state model for the molecular bridge coupled to a single 
optical phonon mode we may write the following expres- 
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FIG. 11: (Color online) Left panel: Schematic of the potential 
energy of the molecular junction versus the occupation num- 
ber no. Right panel: Hysteresis in the current- voltage char- 
acteristics. The solid line corresponds to the transport via 
the unoccupied bridge electronic state, and the dashed line 
corresponds to the transport via the occupied state shifted 
due to the polaron formation. 



sion for the renormalized energy: 

E (n ) = E -^ (32) 

where the electronic population on the bridge no is given 
by: 

no-l/dE f L (E)A L + ME)A R 

W [E — E (no)} + (Al + Ar) 2 

So, as follows from Eq. (|32p the polaron shift depends 
on the bridge occupation n , and the latter is related 
to Eo by the Eq. (|33p . Therefore, the derivation of an 
explicit expression for Eq(tio) is a nontrivial task even 
within the chosen simple model. Nevertheless, it could be 
shown that two local minima emerge in the dependence 
of potential energy of the molecular junction including 
two electrodes linked by the molecular bridge, of the oc- 
cupation number no . These minima are located near 
uq = and no = 1, and they correspond to the neutral 
(unoccupied) and charged (occupied) states, respectively. 
This is illustrated in the Fig. 11 (left panel). These states 
are metastable, and their lifetime could be limited by the 
quantum switching (Mitra 2005, Mozyrsky 2006). When 
the switching time between the two states is longer than 
the characteristic time for the external voltage sweeping, 
one may expect the hysteresis to appear for the states of 
interest live long enough to maintain it. Within the op- 
posite limit the average washes out the hysteresis. When 
the states are especially short-lived this could even result 
in a telegraph noise at finite bias voltage which replaces 
the controlled switching. 

Further we consider long-lived metastable states and 
we concentrate on the I ~ V behavior. Let us for cer- 
tainty assume that the bridge state at zero bias voltage 
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is situated above the Fermi energy of the system and 
remains empty. As the voltage increases, one of the elec- 
trodes chemical potentials crosses the bridge level posi- 
tion, and the current starts to flow through the system. 
The I — V curve reveals a step at the bias voltage value 
corresponding to the crossing of the unoccupied bridge 
level with the energy Eq by the chemical potential. How- 
ever, while the current flows through the bridge, the level 
becomes occupied and, consequently, shifted due to the 
polaron formation. When the bias voltage is reversed the 
current continue to flow through the shifted bridge state 
of the energy E 0} until the recrossing happens. Due 
to the difference in the energies of neutral and charged 
states the step in the I — V curves appears at differ- 
ent values of the voltage, and this is the reason for the 
hysteresis loop to appear as shown in the right panel of 
the figure 11. One could also trace the hysteresis in the 
I—V characteristics starting from the filled (and shifted) 
bridge state. 

Again, the hysteresis loop in the I—V curves may oc- 
cur when both occupied and unoccupied states are rather 
stable which means that the potential barrier separating 
the corresponding minima in the potential energy pro- 
file is high enough, so that quantum switching between 
the states is unlikely. This happens when the bridge is 
weakly coupled to the electrodes (Al : r <C A 2 /fi£l) which 
is an obvious requirement for the involved states (neu- 
tral or charged) to the distinguishable. In other words, 
the broadening of the relevant levels due to the cou- 
pling to the electrodes must be much smaller that the 
polaron shift. As was recently shown (Ryndyk 2008), 
the stronger is the electron-phonon coupling the less be- 
comes the probability of the switching. At large values of 
A the tunneling between the charged and neutral states 
is exponentially suppressed. Also, it was shown that 
the symmetry /asymmetry in the coupling to the elec- 
trodes could significantly affect the hysteresis behavior 
(D'Amico 2008, Ryndyk 2008). For asymmetric junc- 
tions (Ai ^ An) two nearly stable states exist at zero 
bias voltage. When the asymmetry is very strong both 
states could appear stable at one bias voltage polarity 
and unstable when the polarity is reversed. It was sug- 
gested that under appropriate choice of parameters one 
may create a situation when the instability regions for 
the involved states do not overlap. These properties 
give grounds to conjecture that such strongly asymmet- 
ric junctions could reveal memory functionalities, which 
make them potentially useful in fabrication of nanodc- 
vices. 

Among various potentially important properties of 
the electron transport through metal-molecule junctions 
one may separate out the negative differential resistance 
(NDR), that is the decrease in the current I while the 
bias voltage across the molecule increases. The NDR ef- 
fect was originally observed in tunneling semiconducting 
diodes (Kastner 1992). Later, it was viewed in quan- 



tum dots (Grobis 2005, Repp 2005). Several possible 
scenarios were proposed to explain the NDR occurrence 
in the electron transport through molecules. The effect 
could originate from alignment and subsequent disalign- 
ment of the Fermi levels of the electrodes with molecular 
orbitals which happens as the bias voltage varies (Xuc 
1999). Also, the DNR could appear as a Coulomb block- 
ade induced effect (Simonian 2007) and/or due to some 
other reasons. It is likely that different mechanisms could 
play a major part in the NDR appearance in different 
molecular junctions where it was observed so far. 

Here, we discuss the NDR features in the current- 
voltage characteristics which originate from the electron 
coupling with the molecular vibrational modes. As was 
shown (Galpcrin 2005, Ycganch 2007), the polaron for- 
mation could give rise to the NDR. This may happen 
if the polaron shift in the energy of the occupied state 
moves the energy level away from the conduction win- 
dow. As the shift depends on the electronic population, 
and the latter changes as the bias voltage increases, the 
occupied state falls out from the window between the 
chemical potentials of the electrodes at certain value of 
the voltage. If the transport is being conducted via this 
occupied state, then it stops when this voltage value is 
reached. Correspondingly the current value drops to zero 
revealing a distinctive NDR feature. This is a realis- 
tic scenario based on the main features of the electron 
transport through molecules under polaron formation. 

In conclusion, electron-vibrion coupling could lead to 
the formation of a polaron which results in the energy 
differences between occupied and unoccupied states on 
the molecular bridge. At weak interaction between the 
bridge and the electrodes and strong electron-phonon 
coupling these charged and neutral states are metastable 
and could serve for electron conduction. This results in 
such interesting and potentially useful effects as the hys- 
teresis and NDR features in the current-voltage charac- 
teristics. 



VII. Molecular junction conductance and long range 
electron transfer reactions 

Long range electron transfer reactions play an impor- 
tant part in many biological processes such as photosyn- 
thesis and enzime cathalyses (Kuznetsov 1999). Theoret- 
ical and experimental studies of these processes last more 
than four decades but they still remain within a very ac- 
tive research area. In the intramolecular electron transfer 
reactions the electric charge moves from one section of a 
molecule to another section of the same molecule. Long 
range transfer typically occurres in large molecules such 
as proteins and/or DNA, so that these two sections are 
situated far apart from each other. A common setup for 
the transfer reactions includes a donor, a bridge and an 
acceptor. Due to the large distances between the donor 
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(where an electron leaves) and the acceptor (where it 
arrives) typical for the long-range transfer reactions, a 
direct coupling between the two is negligible. Therefore, 
the electron participating in the transfer needs a molec- 
ular bridge providing a set of intermediate states for the 
electron transport. 

Essentially, the intramolecular electron transfer is a 
combination of nuclear environment fluctuations and 
electron tunneling. Electron transfer reactions result 
from the response of a molecule to environmental po- 
larization fluctuations which accompany nuclear fluctu- 
ations. The molecule responds by redistribution of the 
electronic density thus establishing opportunities for the 
charge transfer to occur. The main characteristic of the 
electron transfer processes is the transfer rate K e t which 
is the inversed time of the reaction. The transfer rate for 
the electron tunneling in the molecule interior from the 
donor to the acceptor depends both on the electronic 
transmission amplitudes and the vibrionic coupling. The 
latter provides the energy exchange between the elec- 
tronic and nuclear systems. 

Viewing fast electronic motions in the background of 
the slowly moving nuclei (nonadiabatic electron transfer) 
and applying the Fermi golden rule of perturbation the- 
ory, it was shown that K et could be written in the form 
first suggested by Marcus (Marcus 1965): 

K et = ^\H DA \ 2 (FC). (34) 

Here, the first cof actor is the electron transmission coef- 
ficient, Hda being the effective matrix element between 
the donor and the acceptor. The second term is the den- 
sity of states weighted Franck-Condon factor which de- 
scribes the effect of nuclear motions in the environment. 

There exists a noticeable resemblance between the long 
range electron transfer and the molecular conductance, 
and this resemblance was analyzed in several theoreti- 
cal works (Dahnovsky 2006, Mujica 1994, Nitzan 2001, 
Yeganeh 2007, Zimbovskaya 2002,2003). The fundamen- 
tal part in both molecular conduction and electron trans- 
fer reactions is taken by the electron tunneling inside 
the molecule via the set of intermediate states. In large 
molecules the electron involved in the transfer reaction 
or contributing to the conductance, with a high prob- 
ability follows a few pathways which work as molecular 
bridges. The rest of the molecule do not significantly par- 
ticipate in the transfer/conduction and may be omitted 
from consideration in calculations of the electron trans- 
mission (Jones 2002, Beratan 1992). On these grounds it 
was suggested that the original molecule could be reduced 
to a simpler chain-like structure, and the latter may re- 
place the former in calculations of molecular conductance 
and/or transfer rates (Daizadeh 1997). Such simplifi- 
cation significantly eases computations of the molecular 
conductance and the electron transfer rates in macro- 
molecules. 



The resemblance between the molecular conductance 
and long range electron transfer does not mean that these 
processes are nearly identical. Along with the similarities 
there are substantial differences between the former and 
the latter. For instance, the continuum of states causing 
the charge transport arises from the multitude of elec- 
tronic states on the electrodes in the case of molecular 
conduction, and from the vibrations and fluctuations in 
the environment in the case of electron transfer reactions. 
Also, the driving force which puts electrons in motion 
originates from the external bias voltage applied across 
the junction in the case of molecular conduction. In the 
electron transfer situation this force appears due to the 
electron-vibrionic interactions in the system, or the elec- 
tron transfer reaction starts as a result of photoexcita- 
tion of the donor part of the molecule. The observables, 
namely, the molecular conductance g and/or current I, 
and the transfer rate K et differ as well. 

Notwithstanding these differences the electron trans- 
port through molecular junctions and the electron trans- 
fer could be theoretically analyzed using the common for- 
malism as recently proposed by Yeganeh, Ratner and 
Mujica (Yeganeh 2007). Further we follow their ap- 
proach. The starting point is that one may simulate the 
donor-bridge-acceptor system as some kind of molecular 
junction, shown in the Fig. 12. Here, we represent the 
donor as a single state and this state mimics the 
left lead (assuming the transport from the left to the 
right). The donor is coupled to the bridge and the lat- 
ter is coupled to the continuum of the final states |/) 
simulating the right lead. These couplings are described 
by self-energy terms A^j. For simplicity we assume that 
there exists a single state on the bridge, and this state 
is coupled to a vibrionic mode of the frequency Q. The 
latter mimics environment fluctuations. We assume that 
the "leads" are weakly coupled to the bridge (A^^ <C A) 
which is typical for electron transfer situations. Also, the 
representation of the environment motions by the sin- 
gle vibrionic mode is justified only at low temperatures 
when the thermal energy kT is much smaller than the 
electron-vibrion coupling parameter A. 

Now, we can write the Landauer expression for the 
current flowing through the "junction" using Eqs. ©, 
(|4]). Assuming that the bridge includes only one state we 
simplify Eq. (g]): 

T(E) = 2T r {A 4 GA / Gt} = -2 Ai ^ f Im(G). (35) 

ZAi + A/ 

Here, the subscripts i/f label initial and final states. 
Within the Fermi golden rule regime which allows to in- 
troduce the transfer rate, the bridge must be coupled to 
the final states much stronger than to the initial state 
(A/ > A,) . Therefore: 
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As discussed in Section IV, at weak coupling of the bridge 
to the leads the electron-vibrionic interaction opens the 
set of metastable channels for the electron transport at 
the energies E n = Eq + nhVl (n = 0,1,2...) where Eq 
is the energy of the bridge state with the polaronic shift 
included. The Green's function may be approximated as 
a weighted sum of contributions from these channels: 



E 



P(n) [ej - E Q - nMl + is] 



(42) 



FIG. 12: (Color online) Schematic of the model system used 
in the transfer rate calculations. The initial/final reservoirs 
correspond to the left /right leads in the molecular junction. 



where s — > + , and the coefficients P(n) are probabili- 
ties for the channels to appear given by Eq. fll"8l) . 
Substituting Eq. (g^J) into Eq. (0TJ) we get: 



and the coupling to the final states falls out of the expres- 
sions for the electron transmission and current. Within 
the chosen model this is a physically reasonable result 
for the final states reservoir (the right "electrode" in our 
junction) was merely introduced to impose a continuum 
of states maintaining the transfer process at a steady 
state rate. Also, considering the current flow we may 
suppose that the initial state is always filled [fi (E) = 1] , 
and the final states are empty [ft(E) = 0]. Therefore, 
the current flow through the "junction" accepts the form: 



I=- 



2e 



dEAilm(G). 



(37) 



Both the current and the transfer rate are fluxes closely 
related to each other, namely: K et = I/e. So, we may 
write: 



A" 



2 

irk. 



dEAilm(G). 



(38) 



Now, Aj could be computed using the expression ([7]) 
for the corresponding self-energy term. Keeping in mind 
that the "left electrode" includes a single state with the 
certain energy e$ we obtain: 



A, = Im 



2 S(E 



(39) 



Accordingly, the expression (|38j) for the transfer rate may 
be reduced to the form: 



2 Im[G( ei )] 



(40) 



where represents the coupling between the donor and 
the molecule bridge. It must be stressed that within the 
chosen model, 7j is the only term representing the rel- 
evant state coupling which may be identified with the 
electronic transmission coefficient Hoa in the general 
expression <[34|) for K et . This leaves us with the following 
expression for the Franck-Condon factor: 



(FC) = -~Im[G(ei)]. 

7T 



(41) 



(FC) = P(n)S(AF - nhtt). 



(43) 



71=0 



Here, AF — €i — E n = a — Eo + A 2 / (MY) is the exoergic- 
ity of the transfer reaction, that is the free energy change 
originating from the nuclear displacements accompanied 
by the polarization fluctuations. The effect of the latter is 
inserted via the reorganization term A 2 /(M7) related to 
the polaron formation. The exoergicity in the transfer re- 
action takes on the part similar to that of the bias voltage 
in the electron transport through molecular junctions. It 
gives rise to the electron motion through the molecules. 
In the particular case when the voltage drops between 
the initial state (left electrode) and the molecular bridge 
these two quantities are directly related by |e|V = AF. 

Usually, the long range electron transfer is observed at 
moderately high (room) temperatures, so the low tem- 
perature approximation (|43[) for the Franck-Condon fac- 
tor cannot be employed. However, the expression (|4ip 
remains valid at finite temperatures, only the expression 
for the Green's function must be modified to include the 
thermal effects. It is shown (Yeganeh 2007) that within 
the high-temperature limit (kT > M7) the expression 
for the (FC) may be converted to the well-known form 
first proposed by Marcus: 



(FC) = 



1 



^/AnEpkT 



exp 



(AF 



AE p kT 



(44) 



where E p — A 2 /(Ml) is the reorganization energy. 

While studying the electron transfer reactions in prac- 
tical macromolecules one keeps in mind that both donor 
and acceptor subsystems in the standard donor-bridgc- 
acceptor triad are usually complex structures including 
multiple sites coupled to the bridge. Correspondingly, 
the bridge has a set of entrances and a set of exist which 
an electron can employ. At different values of the tun- 
nel energy different sites of the donor and/or acceptor 
subsystems can give predominant contributions to the 
transfer. Consequently, an electron involved in the trans- 
fer arrives at the bridge and leaves from it via different 



15 



entrances/exits, and it follows different pathways while 
on the bridge. 

Also, nuclear vibrations in the environment could 
strongly affect the electron transmission destroying the 
pathways and providing a transition to completely in- 
coherent sequental hopping mechanism of the electron 
transfer. All this means that proper computation of the 
electron transmission factor Hda for practical macro- 
molecules is a very complicated and nontrivial task. The 
strong resemblance between the electron transfer reac- 
tions and the electron transport through molecules gives 
grounds to believe that studies of molecular conduction 
can provide important information concerning quantum 
dynamics of electron participating in the transfer re- 
actions. One may expect that some intrinsic charac- 
teristics of the intramolecular electron transfer such as 
pathways of tunneling electrons and distinctive features 
of donor/acceptor coupling to the bridge could be ob- 
tained in experiments on the electron transport through 
molecules. For instance, it was recently suggested to 
characterize electron pathways in molecules using the 
inelastic electron tunneling spectroscopy, and other ad- 
vances in this area are to be expected. 

VIII. Concluding remarks 

Presently, the electron transport through molecular- 
scale systems is being intensively studied both theoret- 
ically and experimentally. Largely, unceasing efforts of 
the research community to further advance these studies 
are motivated by important application potentials of sin- 
gle molecules as active elements of various nanodevices 
intended to compliment current silicon based electron- 
ics. Elucidation of physics underlying electron transport 
through molecules is necessary in designing and operating 
molecular-based nanodevices. Elastic mechanisms for the 
electron transport through metal-molecule-metal junc- 
tions are currently understood and described fairly well. 
However, while moving through the molecular bridge, 
an electron is usually affected by the environment, and 
it results in the change of its energy. So, inelastic ef- 
fects appear, and they may bring noticeable changes in 
the electron transport characteristics. Here, we concen- 
trated on inelastic and dissipative effects originating from 
the molecular bridge vibrations and thermally activated 
stochastic fluctuations. For simplicity and also to keep 
this Chapter within a reasonable length we avoided the 
detailed description of computational formalisms com- 
monly used to theoretically analyze electron transport 
through molecular junctions. These formalisms are de- 
scribed elsewhere. We mostly focus on the physics of the 
inelastic effects in the electron transport. Therefore, we 
employ very simple models and techniques. 

We remark that along with nonelasticity originating 
from electron- vibrionic interactions which is the subject 



of the present review there exist inelastic effects of dif- 
ferent kind. For instance, inelastic effects arise due to 
electron-photon interactions. Photoassisted transport 
through molecular junctions was demonstrated and the- 
oretically addressed using several techniques. It is known 
that optical pumping could give rise to the charge flow 
in an unbiased metal-molecule-metal junction, and light 
emission in biased current carrying junctions could occur. 
Also, there is an issue of the molecular geometry. There 
are grounds to conjecture that in some cases the geome- 
try of a molecule included in the junction may change as 
the current flows through the latter for bonds can break 
with enough amount of current. Obviously, this should 
bring a strong inelastic component to the transport, con- 
sequently affecting observables. 

This a common knowledge that electron-electron in- 
teractions may significantly influence molecular conduc- 
tance leading to Coulomb blockade and Kondo effect. 
To properly treat electron transport through molecular 
junctions one must take these interactions into consider- 
ation. Corresponding studies were carried out omitting 
electron-phonon interactions. However, the full treat- 
ment of the problem including both electron-electron 
and electron-phonon interactions is not completed so far. 
There exist other theoretical challenges such as the effect 
of bipolaron formation which originates from the effective 
electron-electron attraction via phonons. 

Finally, practical molecular junctions are complex sys- 
tems, and significant effort is necessary to bring elec- 
tron transport calculations to the result which could be 
successfully compared with the experimental data. For 
this purpose one needs to compute molecular orbitals 
and the voltage distribution inside the junctions, to get 
sufficient information on the vibrionic spectrum of the 
molecule, the electron-phonon coupling strengths and 
electron-electron interactions. One needs a good quanti- 
tative computational scheme for transport calculations 
where all these informations are accounted for. Cur- 
rently, there remain some challenges which have not been 
properly addressed by the theory. Therefore, compari- 
son between theoretical and experimental results on the 
molecular conductance sometimes does not bring satis- 
factory results. However, there are firm grounds to be- 
lieve that further effort of the research community will 
result in detailed understanding of all important aspects 
of the molecular conductance including inelastic and dis- 
sipative effects. Such understanding is paramount to con- 
vert molecular electronics into a viable technology. 
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